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The structure and elastic properties of (5,5) and (10,10) nanotubes, as well as barriers 
for relative rotation of the walls and their relative sliding along the axis in a double-walled 
(5,5)@(10,10) carbon nanotube, are calculated using the density functional method. The 
results of these calculations are the basis for estimating the following physical quantities: 
shear strengths and diffusion coefficients for relative sliding along the axis and rotation of 
the walls, as well as frequencies of relative rotational and translational oscillations of the 
walls. The commensurability-incommensurability phase transition is analyzed. The length 
of the incommensurability defect is estimated on the basis of ah initio calculations. It is 
proposed that (5,5)@(10,10) double-walled carbon nanotube be used as a slider bearing. 
The possibility of experimental verification of the results is discussed. 

I. INTRODUCTION 



In connection with considerable advances made in recent years in the development of nanome- 
chanics, the search for nanoobjects that can be used as mobile elements in nanoelectromechanical 
systems (NEMSs) has become an urgent problem. The possibility of relative motion of the walls 
in multi- walled carbon nanotubes makes it promising to use nanotube walls as mobile elements 
of NEMSs. At present, nanomotors in which the walls of a multi-walled nanotube play the role of 



the shaft and the bush 3|, IJ] and memory cells operating on relative motion of the walls along the 
nanotube axis have been achieved experimentally. A number of NEMSs have been proposed, 
which operate on the re" 

Brownian nanomotor Ml, and a nut-bolt pair 



a varying nanoresistor 



12 



ative motion of carbon nanotube walls, including a gigahertz oscillator [6|], 

9, IQl. Nanoelectromechanical systems, such as 
I. llOl]. a stress nanosensor [11] and an electromechanical nanothermometer 



13l |. have also been proposed. The operating principle of these instruments is based on the 



dependence of the conductivity and the energy of interaction of the walls of a double- walled carbon 
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nanotube (DWNT) on the relative position of the walls on subnanometer scales. Block diagrams 
and operating principles of such NEMSs were 

considered in M. Thus, analysis of the relative 
motion and interaction of carbon nanotube walls is very important for developing and designing 
NEMSs. 

The structure of a nanotube wall is controlled by a pair of integers, chirality indices (n, m) 
corresponding to the coordinates of the carbon lattice vector c = nai + ma2, where ai and a2 are 
the unit vectors of the graphite plane. The segment corresponding to this vector becomes the wall 



circumference when a fragment of the graphite plane is curled to form the nanotube wall 

Neither the interaction between carbon nanotube walls nor the interaction between graphite 
layers have been investigated in details. The results of measurements and theoretical calculations 
of the energy of interaction between graphite layers differ by two orders of magnitude (see 17 1 



and references therein). The shear strengths along the nanotube axis were measured only in a 
few experiments [l], Q]. As regards theoretical studies, the barriers for the relative motion of 



the wal 



s were calculated for a large number of DWNTs using semiempirical 



211], as well as ab initio calculations for certain DWNTs 2^, l23l . \24 . l25l . l26l |. Barriers 



potentials 



for the relative motion of the walls of (5,5)@(10,10) DWNT were determined using all methods 
employed for such calculations. The choice of (5,5)@(10,10) DWNTs was dictated by the following 
considerations: (i) a small number of atoms in the unit cell of a nanotube, which renders it suitable 
for ab initio calculations; (ii) a large height of the barrier for the relative rotation of the walls can 

nn 

be observed only for DWNTs with compatible rotational symmetries of the walls jl9l. I21l|. including 
(5,5)@(10,10) DWNT. 



Certain semiempirical methods 



2ll |. as well as ab initio computational methods [251. l26l| used 
for determining the barriers for the relative motion of the walls of (5,5)@(10,10) DWNT, give val- 
ues of structural parameters and elastic properties of graphite that are in good agreement with 
experiment. Nevertheless, the available data for these barriers in Table 1 show that the results of 
calculations based on different methods differ by several orders of magnitude. Thus, further devel- 
opment of theoretical methods for studying the interaction of layers in carbon nanotubes requires 
that the barriers for the relative motion of the walls be experimentally determined. Comparison of 
experimental and calculated heights of these barriers can be used as a criterion for the correctness 
of the results of calculations. 

Here, we will demonstrate that the results of calculating the elastic properties and structure 
of graphite using the density functional method in the local density approximation with the basis 
consisting of a set of plane waves agree well with experiment. This method is used for calculating 



the elastic properties, the structure, the energy of interaction, and the height of barriers for the 
relative motion of the walls in (5,5)@(10,10) DWNT. The results of these ab initio calculations are 
used to estimate the following physical quantities for (5,5)@(10,10) DWNT: shear strengths and 
diffusion coefficients for the relative sliding along the axis and relative rotation of DWNT walls, 
the frequencies of relative rotational and translational oscillations of the walls, and the energy of 
formation and length of the incommensurability defect. We consider the possible experimental 
study of these quantities. 

It should also be noted that the DWNTs under theoretical investigation in this research were 
obtained using various synthesis methods: in the standard arc for obtaining carbon nanotubes {2?! . 



in the same arc in the presence of hydrogen and a catalyst 

n 



28| . during catalytic decomposition of 



hydrocarbons [29[], and from single- walled nanotu 
to heating jsol or bombardment with electrons 31]. 



3es with a fullerene chain inside and subjected 



II. COMPUTATIONAL TECHNIQUE 



All calculations based on the density functional method in the local density approximation [32 1 
with the basis consisting of a set of plane waves were performed using VASP (Vienna Ab Initio 
Simulation Package) 1331. The interaction of valence electrons with atomic cores was described 

n 

using ultrasoft pseudopotentials [3J]. The maximal kinetic energy of plane waves was 358 eV. Our 
analysis shows that the error in calculating the energy of the system is less than 0.001 meV/atom for 
the chosen cutoff energy for the plane waves. In addition, this error estimate should be insensitive 
to the structure of the system since the basis of plane waves does not contain an error associated 
with superposition of the basis function in contrast to the basis of localized orbitals. 

Summation over the Brillouin zone was performed over special points of the Monkhorst-Pack 
type [35|. The structure and elastic properties of graphite were calculated using a set of 172 
irreducible K points in the Brillouin zone. The barriers for the relative sliding and rotation of the 
walls of (5,5)@(10,10) DWNT were analyzed using a set of nine irreducible K points along the axis 
of nanotubes. 

The structure and elastic properties of graphite were calculated as follows. The lengths of 
vectors a and c of the fundamental translations of the graphite lattice in the plane of the layers 
and in the direction perpendicular to this plane varied in the limits 2.39 A < a < 2.49 A and 
6.57 A < c < 7.90 A, respectively. The calculated dependence of total energy Ug of the system 
on the length of the fundamental translation vectors of the graphite lattice was approximated by 



a third-degree polynomial, 

Table 1 Barriers ISUz and ISU^f, for relative sliding of the walls along the axis and relative 
rotation of the walls about the (5,5)@(10,10) DWNT axis per atom of the outer wall and ratio 
7f, = AU^/AUz of these barriers 



Reference 


AUz, meV/atom 


AC/^, meV/atom 


76 




0.008 


0.025 


3.13 




0.00745 


0.0144 


2.90 


[20]= 


7.5 


8.7 


1.16 


[24]"^ 


0.128 


0.438 


3.47 


[23]-^ 




1.2 




[22]*^ 


0.35 


0.78 


2.26 


[25]^ 


0.125 


0.259 


2.08 


[26]5 


0.138 


0.407 


2.85 


Our data= 


0.005±0.001 


0.205±0.002 


41 



Note: " Lennard-Jones potential, optimized structure of the walls; Lennard-Jones potential, nonoptimized 
structure of the walls; Crespi-Kolmogorov potential, optimized structure of the walls; tight binding 
method; density functional method, local density approximation, basis set is formed by plane waves; 
density functional method, local density approximation, basis set is formed by pseudow ave pdpp functions; ^ 
density functional method, local density approximation, basis set is formed by pseudowave dddd functions. 

Table 2 Coefficients of approximation of the dependence of the total energy of graphite on the lengths 
of fundamental translation vectors of the lattice by third-degree polynomial ([T]) 



Coefficient 


Value, eV/unit cell 


Coefficient 


Value, eV/unit cell 


Don 


-40.6169±0.00004 


D20 


-49.0±0.5 


Dsn 


0.0580±0.0003 


D02 


0.036±0.003 


Do3 


46.98±0.02 


Dii 


0.011±0.033 


D21 


-0.069±0.003 


D12 


-0.0223±0.0002 



Ug = Doo + D3o{c - cof + Do3{a - gq)^ + D2i{c - co)^{a - ao) + £'12(0 - co)(a - aof + 

+ D2a{c ~ co)^ + £'02(0 - ao)^ + £'ii(c - co)(a - oq), (1) 

where oq and cq are the equilibrium lengths of the fundamental translation vectors of the graphite lattice. 
Table 2 gives the values of the coefficients obtained as a result of approximation, while Table 3 gives the 
values of equilibrium lengths of the fundamental translation vectors of the graphite cell. The elastic constants 
of graphite were obtained using expression ([T|) for the total energy of the system, 



where Vq is the equiUbrium volume of the unit ceU of graphite. Table 3 contains the elastic constants of 
graphite. In the same table, our results of calculating the structural and elastic properties of graphite with 
the density functional method in the local density approximation are compared with the results obtained in 
[3^ and the experimental data described in 



This comparison shows that the accuracy of our calculations 
corresponds to the typical accuracy of calculations by the density functional method. 

Thus, our calculation method gives values of the lattice parameters and elastic properties of graphite 
close to experimental ones. We believe that this method can also be used for studying the interaction and 
relative motion of carbon nanotube walls. 



Table 3 Equilibrium lengths of the fundamental translation vectors of the lattice and elastic constants 
of graphite 





Our result 


[36] 


Experiment (300K) 


Oo, A 


2.44140±0.00001 


2.440 


2.4608±0.0016 


Co, A 


6.586±0.002 






co/oo 


2.698 


2.74 


2.725±0.001 


Cii +C12, GPa 


1319 


1283 


1240±40 


C33, GPa 


23.7 


29 


36.5±1 



Table 4 Structural parameters and elastic properties of the (5,5) and (10,10) walls: R is the radius, tg 
is the equilibrium length of a unit cell, and Es is the Young modulus 



Wall 


(5,5) 


(10,10) 




our result 


available data 


our result 


available data 


R, A 


3.40 


3.408[26] 


6.77 


6.748 [26] 


ts, A 


2.439 


2.4434[26] 


2.443 


2.4420f26] 


Es, TPa 


1.24 


0.95[38], 0.96[39], 


1.31 


0.92[38], 0.98[40], 






1.03f26], 1.06f37] 




1.10[26], 1.24[41] 



III. RESULTS AND DISCUSSION 

First, we calculated the equilibrium structure of isolated (5,5) and (10,10) walls. The structure of each 
wall was optimized for various lengths t of the unit cell of the wall. The dependence of the total energy Ut 
of the unit cell on its length was interpolated by Hooke law, 



(3) 



where Us and tg are the total energy and length of the unit cell, which correspond to the equilibrium 
structure, and is the elastic coefficient of a wall of length tg. The values of Us, tg, and /i were calculated 
using the least squares method. The Young modulus of the wall was determined in most publications using 
the expression 

= it- <"> 

where R is the radius of the wall corresponding to the equilibrium structure and w is the effective thickness 
of the waU: 3.4 A bfi 



37 



38 



38 



4C 



4l|. Table 4 gives the calculated values of the structural parameters 
and the Young modulus of the (5,5) and (10,10) walls, as well as the corresponding data from the literature. 
It should be noted that the calculated lengths of the bonds are in conformity with the results of other 



calculations based on the density functional method 
match those calculated using the density functional 



4C 



4l|. 



26l ■ while the calculated values of the Young modulus 



39| , Hartree-Fock [37i| , and tight binding methods 



The difference in the unit cell lengths of the (5,5) and (10,10) walls, Isi — t2 — ti — 0.004 A, is three 
orders of magnitude smaller than the unit cell lengths. However, a wall with a smaller unit cell length can 
be extended, while that with a larger unit cell length can contract as a result of the interaction between the 
walls H]. Here, the energy of interaction between the walls is calculated for DWNTs with walls having an 
optimized structure and with a unit cell of the same average length ta — {t2 + ii)/2. 

We determined the energy of interaction between the walls as the difference between the total energy of 
the DWNT and the total energies of isolated walls. To analyze the characteristics of the relative motion of 
nanotube walls, we must calculate the dependence of energy U of interaction between two adjacent walls 
on the coordinates describing the relative position of the walls (angle (f) of the relative rotation of the walls 
about the nanotube axis and length z of the relative displacement of the walls along this axis) . 

In accordance with the symmetry of (ni, 7ii)@(ri2, 712) DWNTs, the unit cell of the potential re- 
lief of the interwall interaction energy of such DWNTs is a rectangle with sides 5z = ta/2 and 6^ = 
T:GCD{ni, /nin2, where GCD{ni,n2) is the greatest common division of numbers rii and n2 42]. For 
(5,5)@(10,10) DWNT, we have 5^ = tt/IO. The interwaU interaction energy U{z, (f>) for (5,5)0(10,10) DWNT 
was calculated for 100 relative positions of the walls, corresponding to a unit cell of the potential relief. We 
found that potential energy Uint of the interaction between layers, which corresponds to the minimum of 
function U{z, 0), is 14.5 meV per atom of the outer wall. This result is in good agreement with the exper- 
imental values of the energy of interaction between graphite layers (35±10 meV/atom) tl7.j and between 
the walls of multi- walled nanotubes with large diameter (20-33 meV/atom) [2], as well as with the results 
of calculation by the density functional method in the local density approximation with the basis set of 
pseudowave dddd functions (23.83 meV/atom of the outer waU of a (5,5)@(10,10) DWNT) 3- 

The relative positions of the walls corresponding to a higher symmetry of a DWNT are singular points of 
interwall interaction energy U{z, (f) (extrema or saddle points) For such relative positions of the walls, 
some of the U2 axes (binary axes) coincide. The U2 axes are perpendicular to the principal DWNT axis and 



pass through the middles of all bonds and through the centers of all hexagons of the structure of the wall. 
According to our calculations, the U2 axes passing through the middles of tilted bonds of the (5,5) wall 
and the middles of the bonds lying in the plane perpendicular to the DWNT axis of the (10,10) wall layer 
coincide for the relative position of the walls corresponding to the minimum of the energy of interaction 
between the walls. The same relative position of the walls corresponding to the minimum of energy[/(z, 0) 



was obtained using semiempirical calculations [43|. The relative position of the walls corresponding to the 
energy minimum in our calculations is shifted along the z axis by a half-period Sz relative to the position 
obtained from calculations based on the de nsity functional method in the local density approximation with 
the basis set of pseudowave dddd functions 26 1. 



Table 1 (see above) contains the calculated values of barriers AUz and AC/^ for relative sliding of the walls 
along the DWNT axis and relative rotation of the walls (between their relative positions corresponding to the 
interaction energy minimum), respectively, as well as ratio 7b = AU^/ AUz of these barriers and analogous 
results from the literature obtained with the help of semiempirical methods and ab initio calculations. 

In accordance with calculations based on the density functional method {2^ and semiempirical calcu- 
lations 2l|, the energy of interaction between commensurate nonchiral DWNT walls (ni, rti)@(n2, ^2) and 



(ni, 0)@(n2, 0) can be interpolated by the following expression: 

TKA ^ TT ^^t' /^27r \ AUz (2n \ 

U{(t>,z) = UQ 2^°^\6^) 2~^°^\S~^)' ^' 

where Uq is the mean energy of interwall interaction. Henceforth, we will normalize the values of Uq, AUz 
and AU^ to an atom of the outer wall. In accordance with our calculations, the shape of barriers for relative 
sliding and rotation of the walls can also be interpolated by cosinusoids to within the computation error; 
however, the heights of the barriers for sliding or rotation depend on the relative rotation or shear of walls, 
respectively, along the DWNT axis. Table 1 contains our results for barrier heights for relative rotations and 
sliding of the walls at the instant of their passage through the relative position of the walls corresponding 
to the minimum of the interwall interaction energy. However, the barrier height at the instant of passage 
through the relative position of the walls corresponding to the energy maximum of interwall interaction is 
2% lower as compared to that for the passage through the minimum, and the barrier height for the sliding 
is lower than the computational error (0.001 meV/atom). 

The high value of the barrier ratio for relative rotations and sliding of the walls (see Table 1) also indicates 
a high ratio of threshold forces for relative rotations and sliding of the walls. Consequently, we believe that 
this nanotube can be used as a slider bearing. 

A telescopic extension of the inner wall gives rise of a capillary force Fc pulling the inner wall back into 
the outer wall. The average value of this force is defined as 

- ) 

where 4n2 is the number of atoms in the unit cell of the outer wall. The average capillary force for interwall 



interaction energy values Uq and length ta of the unit cell of the nanotubes calculated here is {Fc) — 0.380 
nN. 

Our values of barriers AUz and AU^ for the relative motion of DWNT walls can be used for calculating 
the values of a number of physical quantities characterizing the interaction and relative motion of DWNT 
walls. Expression ^ for the dependence of interwall interaction energy U (0, z) on the relative position of 
the walls defines the threshold static friction forces Fz and F^ for the relative sliding of the walls along the 
nanotube axis and their relative rotation, respectively: 

_ 4n2TrLo^AUz _ 4n2'!rLo^AUz , , 

where Lqv is the overlap length of the walls and Rm is the radius of the movable wall. The shear strengths 
for relative sliding of the walls along the nanotube axis and their relative rotation are defined, respectively, 
as 

Mz = ^, = ^, (8) 

where S is the area of the overlap surface of the walls, 

S = 2ttLo. = 7rLo.{Ri + R2), (9) 

i?i and i?2 are the radii of the inner and outer wall, respectively. In accordance with our calculations, the 
values of shear strengths for (5,5)@(1G,10) DWNT are Mz = 1.0 ± 0.2 MPa and Af^ 51.5 ± 0.5 MPa, 
respectively, for the relative sliding of the walls along the nanotube axis and their relative rotation. In 
calculating the shear strengths for the relative rotation of the walls, we assumed that the inner wall is 
mobile. 

The equality of forces Fc — Fz acting on the movable wall controls the minimal overlap length i,„ of the 
walls for which the static friction force prevents the inner wall from being pulled by the capillary force: 

According to our calculations, = 110 ± 20 nm for (5,5)@(10,10) DWNT. 

To find the frequency of small rotational oscillations and relative vibrations of the walls along the DWNT 
axis, we approximate expression for the dependence of the interwall interaction energy on their relative 
position by a parabolic potential well in the vicinity of the minimum Um of this dependence. 



U{Z') = [/„ + U{C^') = Urn + (11) 

where z' and 0' are the coordinate of displacement relative to the energy minimum of interwall interaction. 
For (m, ni)@(n2, n2) DWNTs, we have 

_ 8n2LoyAUz ( 27r^^ _ 8n2LoyAU^ f nirt2 \^ 

\TJ ' \GCD{n,^n2)) ' ^^^^ 



Frequencies v^p and Vz of small relative rotational and axial (along the DWNT axis) oscillations of the 
walls are defined as 



respectively, where M and J are the reduced mass and moment of inertia of the walls. 

In the case when only the inner wall of a DWNT is mobile, its moment of inertia is Ji = MiRf, 
where Mi = AnimoLi/ta is the mass of the inner wall, Li is its length, and toq is the mass of a carbon 
atom. Analogously, if only the outer wall is mobile, its moment of inertia and mass are defined as J2 = 
M2R2 and AI2 = 4:n2moL2/ta, where L2 is the length of the outer wall. In the case when both walls are 
mobile, the reduced moment of inertia and the mass of the system are given by J12 = JiJ2/{Ji + J2) and 
A/12 — A/iM2/(Mi + Af2). The latter case can be encountered, for example, in zero-gravity conditions in a 
space-based laboratory. 

The frequencies of rotational and axial relative oscillations of the walls are independent of the DWNT 
length when the lengths of both walls are identical, or when a shorter wall is mobile and the longer wall is 
fixed. In these cases, the overlap length of the walls coincides with the length of the mobile wall, and the 
expressions of the frequencies of relative axial oscillations of the walls assume the form 



1 2712 AUz 1 2AUz 

1^z,l^—\ , Vz.2 = —\ , (14) 

ta V "iTlo ta V "IQ 



1 2{m+n2)AUz 

ta y niTOo 

Here, i^z,i, Vz.2i and Vz,i2 correspond to the inner, outer, and both mobile walls. The expressions for 
frequencies i^^.i, ^^0,2, and 1^0,12 of relative rotational oscillations of the walls corresponding to the inner, 
outer, and both mobile walls, respectively, assume the form 



n2 nin2AU^ _ nin2 AU^ 

"''''^ ~ TrGCD{nin2)Ri\l 2mo ' ^''''^ ~ nGCD{nin2)R2 \ 2mo' ^ ' 



n2 ni{niRl + n2Rj)AU^ 

"'''■'^^ ~ TrGCD{nin2)RiR2 V 2mo ' ^ ^ 

As a result of our ab initio calculations, we obtained the following values of frequencies of relative 
oscillations of the walls: Vz^i = 1.74 ± 0.17, i^z,2 = 1-22 ± 0.12, iyz,i2 = 2.1 ± 0.2 = 12.66 ± 0.06, 
z^0,2 = 4.50 ± 0.02 and 1^^,12 = 13.44 ± 0.07 cm^^. 

If the energy of thermal motion of a short mobile wall of length L is much lower than the height of 
the barriers for the sliding of this wall along the DWNT axis {kT ^ An2AUzL/ta) and/or for its rotation 
[kT <^ An2AU^L /to), then hopping diffusion of the shorter wall along the DWNT axis and/or its hopping 
rotational diffusion, respectively, takes place Q, [l^l- In this case, diffusion coefficients Dz and for the 
motion of the mobile wall along the axis and for rotational diffusion are defined as 



1_ 2 f ^n2MJ,L 



i^. = 2"A^exp^-^^j (18) 

respectively, where and are the preexponential factors in the Arrhenius formula for the frequency of 
the jump of the mobile wall between equivalent minima of interwall interaction energy. 

Simulation based on the molecular dynamics method shows that the preexponential factor in the Ar- 
rhenius formula for the frequency of the jump of a C6o@C24o nanoparticle shell between equivalent energy 



minima of interwall interaction is = 540 ± 180 GHz 



4J] . The frequency of small relative rotations near the 



minima of interwall interaction energy calculated here for such a nanoparticle is = 60 GHz for the same 
shape of the shells and the same potential of interaction between atoms of the shells. Thus, ratio Vl/v oi the 
preexponential factor in the Arrhenius formula to the frequency of small oscillations is approximately 10. We 
believe that ratio ^l/v for other carbon nanostructures with embedded graphite layers (in particular, carbon 
nanotubes) is on the same order of magnitude. We used the value oi ^/v ^ IQ for estimating the diffusion 
coefRcients of walls in (5,5)@(10,10) DWNT on the basis of our ab initio calculations. The values for the 
diffusion coefficients of a short mobile wall of length L ~ 100 nm along the DWNT axis at a temperature of 
300 K are Dz.i ~ 8 ■ 10^^^ ni^/s and ^ ~ 6 ■ 10^^^ m^/s for the inner and outer wall, respectively. 

Under the action of force Fd directed along the DWNT axis, the shorter wall may drift along the longer 
one at a velocity v — B^Fd, where is the mobility of the shorter wall for its motion along the DWNT axis. 
In the case of a small force {FdSz "C 2AU, where AU is the barrier height for the motion of the shorter wall 
along the DWNT axis), mobility Bz is related to the diffusion coefficient by the Einstein relation Dz = kTBz 



10|. According to our estimates, for the drift of a wall 100 nm in length at a velocity of 0.1 m/s at a 
temperature of T = 300 K, forces of Fd,i ~ 50 pN and Fd,2 ~ 70 pN must be applied to the inner and outer 
walls, respectively. 

Rotational diffusion for walls of length L = 100 nm at a temperature of 300 K in (5,5)@(10,10) DWNT 
is ruled out in view of the large height of the barrier for the relative rotation of the walls. For a mobile wall 
of length L = 10 nm, the values of rotational diffusion coefficients at 300 K are ZJ^.i « 2 • 10'' rad^/s and 
-D0_2 ~ 6 • lO'^ rad^/s for diffusion of the inner and outer wall, respectively. 

A theory of the commensurability-incommensurability phase transition in DWNTs with nearly commen- 
surate walls has been developed recently [llj] . In accordance with this theory, such a DWNT can be either in 
a commensurate phase with identical lengths of the unit cells of the walls, or in an incommensurate phase, in 
which long segments with practically commensurate walls alternate with relatively short incommensurability 
defects. The phase of the DWNT is controlled by the value of the incommensurability parameter 

2n2AUz MiM2 .„„s 

Mi2(Ai)^ 



where fix and ^2 are the elastic coefficients of the inner and outer wall, respectively, with the length of a 
DWNT unit cell. 

It should be noted that the incommensurability parameter is the ratio of the difference in energy of 
interwall interaction for the commensurate phase and a completely incommensurate state to the elastic 
energy of the DWNT in the commensurate phase (both values of energy are normalized to a DWNT unit 
cell). The DWNT is in a commensurate phase ii h > he = and in an incommensurate phase if h < he- 
In accordance with our calculations, h = 0.13 ± 0.03 for (5,5)@(10,10) DWNT. Thus, this DWNT is in 
an incommensurate phase. In this phase, the DWNT has the following structure: commensurate regions 
with the same length of the unit cells of the walls alternate with incommensurability defects. In each 
incommensurability defect, the number of unit cells in the wall with a shorter equilibrium length of the 
unit cell is greater by unity than the corresponding number in the other wall (incommensurability defects 
are analogous to dislocations in a crystal). The expression for the length of an incommensurability defect 
assumes the form [ll| 

Our calculations give an estimate of Id ~ 146 ± 21 nm. 



IV. CONCLUSIONS 



Let us consider how to experimentally verify our results. At present, the upper boundary of the shear 
strength for the relative sliding of the walls relative to the nanotube axis determined with the help of atomic 
force microscopy is Mz < 0.04 In most cases, the barriers for the sliding of adjacent walls (and, hence, 
the corresponding shear strength) are negligibly small (due to the incompatibility of translational symmetry 



19| and to incompatibility of helical symmetries of the walls for 



of the walls for incommensurate walls 

commensurate chiral walls [l^H, [3]). In experiment the shear strength was measured only for a single 
pair of walls of a multi-walled nanotube, which was characterized by a smaller value of this strength as 
compared to other pairs of adjacent walls. The chirality indices of the walls were not determined. Thus, 
we believe that the experimental value of the upper boundary of the shear strength for the relative sliding 
of the walls corresponds to incommensurate or commensurate chiral walls. However, the barrier for the 
relative sliding of adjacent walls may attain an a ppr eciable value only for commensurate nonchiral walls (in 



particular, for (5,5)@(10,10) DWNT) 3, Q, [2O , l2l| . The shear strengths calculated here for (5,5)@(10,10) 



DWNT exceed by several orders of magnitude the corresponding boundary of this value (0.04 MPa) measured 
by an atomic force microscope for the relative sliding of adjacent walls with indeterminate chirality indices. 
Consequently, the shear strengths for this DWNT can easily be determined experimentally. In our opinion, 
the minimal overlap length of the walls in the case of the telescopic extension of the inner wall, for which the 
static friction force prevents the inner wall from being pulled by the capillary force, can also be measured 
using atomic force microscopy. 



The DWNT images obtained by transmission electron microscopy demonstrate that the interwall interac- 
tion may result in extension of one wall and compression of the other wall accompanied by the formation of 
commensurate regions j3| ■ We believe that transmission electron microscopy can also be used for observing 
incommensurability defects in DWNTs. 

The calculated frequencies of relative oscillations of the walls can be measured using terahertz spec- 
troscopy or Raman spectroscopy. It should be noted that Raman spectra have recently been obtained for 
isolated nanotubes 46]. In such measurements, the frequencies of relative oscillations of the walls should 
be distinguished from the frequencies of other DWNT modes corresponding to those in isolated walls. In 
accordance with our calculations, the modes of (5,5)@(10,10) DWNT with the lowest fre quen cies, which 
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Thus, the 



correspond to the modes existing in isolated walls, have frequencies exceeding 40 cm"-'^ 
frequencies of relative oscillations of the walls in (5,5)@(10,10) DWNT correspond to another frequency 
range as compared to the frequencies of the remaining modes of this DWNT. 

Since the barriers for the relative motion of the walls are controlled by the chirality indices of the walls 
0, H, H, 0, H, [25 , 26 1, the measurements of the shear strengths for the relative motion of the walls and 



the frequencies of their relative oscillations should be accompanied by determination of the chirality indices 

of the walls. It was shown that the chirality indices of both walls of a DWNT can be determined using 

f 

electron diffractomctry 49]. We believe that the characteristics of (5,5)@(10,10) DWNT obtained here on 
the basis of ab initio calculations of the interwall interaction energy can be determined by contemporary 
experimental techniques. Such measurements would facilitate progress in developing theoretical methods 
for studying the interwall interaction in nanotubes, as well as in developing NEMSs based on the relative 
motion of the walls of nanotubes. 
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